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We theoretically examine three-well interferometry in Bose-Einstein condensates using adiabatic 
passage. Specifically, we demonstrate that a fractional coherent transport adiabatic passage protocol 
enables stable spatial splitting in the presence of nonlinear interactions. A reversal of this protocol 
produces a coherent recombination of the BEC with a phase-dependent population of the three 
wells. The effect of nonlinear interactions on the interferometric measurement is quantified and 
found to lead to an enhancement in sensitivity for moderate interaction strengths. 
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I. INTRODUCTION 

Since the creation of the first Bose-Einstein conden- 
sate (BEC) [1, 2 , cold atom experiments have provided 
a powerful platform for the study of macroscopic quan- 
tum states [31 H] , emulation of solid state physics [5^, W , 
and insight into many-body quantum phenomena. The 
utility of BECs stems from the experimental freedom to 
control many system parameters, including the interac- 
tion strength, degrees of freedom, size, and shape of the 
BEC. 

One enticing proposal is to use BECs for interferome- 
try [THQ]. BEC interferometers would have many advan- 
tages over their optical counterparts. Trapped-atom in- 
terferometers can be sensitive to changes in mass, charge, 
magnetic moment and polarisability. As with optical in- 
terferometry, a BEC interferometer involves the spatial 
splitting, followed by the generation of a relative phase 
difference between split components and then coherent 
recombination of the quantum state. Performing these 
operations in a way that is relatively insensitive to small 
errors in the implementation is a necessary precondition 
for a reliable interferometer. 

The conventional two-well splitting of a BEC is highly 
sensitive to atom- atom interactions, where phase diffu- 
sion and Josephson oscillations lead to a loss of phase 
resolution p!Qm2 . Beam splitting via laser-induced adi- 
abatic passage has the advantage of being robust to 
changes in experimental parameters. As first proposed 
by Marte et at [13 a variation of the efficient and ro- 
bust three-state adiabatic process known as stimulated 
Raman adiabatic passage (STIRAP) [MHIZ], called frac- 
tional STIRAP (f-STIRAP) [M [19], can be used to 
generate any preselected coherent superposition of two 
atomic states, |1) and |3), via an intermediate excited 
state, 1 2). Electromagnetic pulses are used to couple 
states |1) to 1 2) and |2) to |3), characterised by coupling 
parameters Vti2 and (^23- As in STIRAP, (^23 precedes 
but unlike STIRAP where 1^23 vanishes first, here 
the two pulses vanish simultaneously while maintaining 
a constant ratio of amplitudes. The ratio of probabil- 
ity amplitudes of the resulting coherent superposition of 



states |1) and |3) is proportional to the ratio ^^23/^12- 
Hence stopping at ^^23/^12 = 1 can create an atomic 
beam splitter, as demonstrated experimentally by Weitz 
et al [ini. 

Here we propose an alternative method for the spa- 
tial splitting and coherent recombination of a BEC 
based on the ideas underpinning f-STIRAP. In three- 
well atomic [21 and electronic quantum dot [22 sys- 
tems, the coherent spatial transport of single particle 
quantum states is known as coherent tunnelling adiabatic 
passage (CTAP). Recent work has shown that this princi- 
ple can be extended to interacting many-body quantum 
systems such as BECs, both in the quantum [23 and 
semi-classical mean-field j24H3Ql limits. 

Analogously to f-STIRAP, fractional CTAP (f-CTAP) 
allows for the creation of a coherent spatial superposition 
of eigenstates of wells 1 and 3. For photons, f-CTAP has 
been demonstrated experimentally in three-channel op- 
tical waveguides [31 . Consider the initial state with the 
BEC confined to well 1, schematically shown in Fig.jlja). 
The atomic population of the BEC is transported from 
well 1 to 3 through adiabatic changes to the tunnelling 
rates between the wells, and is halted once the BEC is 
equally split [see Fig.jljb)]. At this point, one component 
of the BEC can be made to interact with some system of 
interest for time, r. In the case where a phase difference 
accumulates between the two states of the superposition, 
this phase can be determined via an interferometric mea- 
surement, as demonstrated in two- well systems [Q] [32] 
by releasing the BEC from the three-well system. In 
this work, f-CTAP is proposed to be also used in the re- 
combination stage of the interferometer as reversing the 
f-CTAP protocol results in a phase-dependent popula- 
tion of the three wells [see Fig. [ijc)]. Interferometric 
f-CTAP needs to be performed on a timescale that is 
much longer than non-adiabatic methods. However, the 
tunnelling interaction between the split BECs maintains 
mutual coherence throughout the splitting and recombi- 
nation processes. As with CTAP, the f-CTAP protocol 
has the advantage of being robust to variations in trap 
parameters and pulse time. 

To explore the dynamics of our three-well system, we 
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FIG. 1: Schematic representation of three- well system: at (a) 
t = 0, at (b) tp/2 < t < tp/2 + r where r is the hold time 
required for phase accumulation, and at (c) t = tp + r. The 
system consists of two parallel, repulsive Gaussian barriers 
embedded in an ambient harmonic trap, dividing the system 
into three wells. At t = - (a) - the BEC resides in well . 
At t = tp/2 - (b) - the BEC is split into two equal compo- 
nents residing in wells 1 and 3. To perform the interferometric 
measurement, the tunnelling rates are then kept constant for 
some hold time, r, during which a relative phase difference 
may be accumulated. At t = tp + r - (c) - the two BECs are 
recombined which leads to a phase-dependent population in 
well 1. 



interaction f-CTAP model can be directly mapped onto 
a corresponding Gross-Pitaevskii equation (GPE) [26 1 133 1 
[34], however the key effects can all be studied using the 
three-mode model without the added complexity of a full 
GPE simulation. 



II. NONINTERACTING MODAL 
APPROXIMATION 

Consider a three-well system. For sufficiently large 
wells, the system is accurately described by the modal 
approximation [TTl ES] [26l [28] . Each well is approximated 
by a single mode, giving: 



(1) 



where the amplitude of each mode j is expressed as t/jj = 
yJW^e^^^ ^ with and 0j being the modal occupation 
and phase respectively. The system is normalised such 
that Xlj=i^i(^) — ^- ^^is basis, the Hamiltonian of 
the three-well system is: 



H = hn 




(2) 



where Qij is the dimensionless tunnelling rate between 
modes i and j, determined by the wavefunction overlap. 
The dimensionless on-site interaction energy per particle 
is Uj = Ej -\- gNj . Ej is the groundstate energy of well 
j, in units of hQ^ with Q being the maximum tunnelling 
rate and g the dimensionless parameter describing the 
nonlinear atomic interactions. 

We first describe the f-CTAP protocol for BEC split- 
ting and recombination in the noninteracting regime, 
^ = 0. 



A. BEC Splitting with Fractional CTAP 



employ two approximations of varying levels of sophisti- 
cation. The first is a three-mode approximation without 
interactions pSJ [26l [28] . In this limit we show analyti- 
cally that a BEC initially residing in well 1 can be split 
between wells 1 and 3, and recombined to give a phase de- 
pendent density in well 1 [see Fig.[l]. To understand the 
role of interactions we use a nonlinear three-mode treat- 
ment. The presence of interactions gives rise to a window 
where canonical CTAP can occur and also a regime where 
more sensitive interferometric sensing can be observed. 
By solving the corresponding classical equations of mo- 
tion, we identify the nonlinear eigenstates of the system, 
identifying the bifurcation point where extra nonlinear 
states appear, and investigate the stability of the CTAP 
state. 

The nonlinear dynamics that are obtained from the 



In the absence of nonlinear atomic interactions = 0) 
and for Ef = ^3 = 0? £'2 = ^5 eigenstates of the 
Hamiltonian, Eq. (??), are: 



^12^3 



\J 1^23 + ^12 



(3) 



V7(7±A) 



721^12^1 



A±7 



^2 + V^r^23^3 



(4) 



where 7 = ^J^^^^{Vfi^ 
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The 1^0 eigenstate 
has no population in the centre well for all A and in the 
limit of ^^12/^23 ^ 1 the atomic population is confined 
entirely to well 1. 

Consider the scenario where the BEC initially resides 
in well 1 (^^12/^23 = 0). Adiabatic transport through the 
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FIG. 2: Ideal three- well system with A = and ^ = 0: (a) 
Proposed pulsing scheme: Eqs. (??) and (??). (b) Energies of 
the eigenmodes: Do and D±. (c) Evolution of the occupation 
of the three wells: N\ (solid curve) and N-^ (dashed curve), 
(d) Adiabaticity parameter, cd^. In these figures we have 
assumed that r = 0. 



1^0 state such that the BEC is equally split between wells 
1 and 3, at t = requires a smooth transformation of 
Vt\2 = to 1^12 = 1^23 = 1/2. A pulsing sequence, for 
r = that meets this criterion is: 



^23 (^) = - + - COS [TTt/tp] , 



(5) 
(6) 



as shown in Fig. |2|a). The mode energies for this pulse 
sequence are shown in Fig. [2|b) with the population of 
the three wells shown in Fig7|2|c) (where adiabatic puls- 
ing has been assumed). At N\ = A/3 =0.5 which is 
the condition we require for sensing. 

The robustness of the splitting protocol relies on adi- 
abatic changes to the system tunnelling. For noninter- 
acting BECs, the adiabaticity is quantified by the adia- 
baticity parameter: 



max 



{Do;t\H{t)\D,;t) 



\Eo{t) - E,{t)\' 



(7) 



which is shown in Fig. |2|d) through the protocol. Adia- 
batic transfer requires tp to be chosen such that a <C 1 
over the entirety of the protocol. 



B. Phase Interferometry with Fractional CTAP 

Once the BEC is split, interaction with a target sys- 
tem can induce a relative phase difference between the 



two components. Since 11 12 = ^23 = 1/2 at t = tp/2, 
phase information is shared between the BECs in wells 
1 and 3. This enables the mutual coherence between the 
two BECs during splitting to be maintained, inhibiting 
phase diffusion. To enable the desired phase accumula- 
tion, tunnelling between wells 1 and 3 needs to be sup- 
pressed. This can be achieved via a variety of means, 
such as (i) applying a blocking laser field at the middle 
point of well 2 or (ii) reducing 1^12 and 1^23 adiabatically 
to zero. After the decoupling of the BECs, phase diffu- 
sion will occur. The BECs can be recoupled some time 
later by, for example, (i) switching off the blocking laser 
field or (ii) adiabatically increasing 1^12 and 1^23 to 1/2. 
We assume that the process of decoupling and recoupling 
occurs over timescale r. Considering the case where the 
two BECs accumulate a phase difference (j) over r, the 
state of the system is then: 



^5^/2 + r,( 
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(^1 - e'^^s 



(8) 



The acquired phase difference during r moves the sys- 
tem out of the 1^0 state, and is now in a superposition 
of all three eigenstates. Hence, the resulting change in 
eigenstate populations for a given (j) is: 



N{D^) = sin^((^/2)sin^(9, 
N{D_) = sin^(0/2)cos^ l9, 
N{Do) = cos\^/2), 



where 



arctan 



2/A 



(9) 
(10) 

(11) 
(12) 



Adiabatic evolution of the system through to t = tp -\- r 
returns the barriers back to their initial configuration, 
where 1^12 = and 1^23 = 1. At the end of the splitting 
and phase accumulation (t = + r) , the I^o eigenstate 
is simply the state in well 1. Thus the final population 
in well 1 is given by: 



Ni{tp- 



N{Do) = cos\^/2). 



(13) 



Thus, the phase difference accumulated during r man- 
ifests as a population imbalance at the end of the pro- 
tocol. Density measurements on the final state of BEC 
then allow determination of the phase, (j). 



III. NONLINEAR THREE-MODE 
APPROXIMATION 

Controllable nonlinearities are one of the defining fea- 
tures of BEC dynamics as compared to analogous pho- 
tonics systems. Accordingly we now consider the effect 
of a non-zero g on the f-CTAP interferometric protocol. 
Inter-particle interaction in BECs gives rise to a raft of 
interesting phenomena, for example quantum phase tran- 
sitions [6 , solitons |35H38] . and entanglement generation. 
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However, in the context of BEC interferometry, these 
interactions are in general problematic, as they induce 
phase diffusion. These interactions manifest as a non- 
linearity in the mean- field description, leading to highly 
non-trivial dynamics. We analyse the behavior of the 
system by mapping the nonlinear Schrodinger equation 
to a classical Hamiltonian. The nature of the stationary 
states in the classical model determines the behaviour of 
the protocol in the adiabatic limit. This analysis is com- 
plimented by numerical simulation. Simulation of the 
protocol confirms the conclusions of the classical analy- 
sis in the large tp limit, and reveals several interesting 
phenomena for finite tp. Interestingly, we find that the 
presence of nonlinearities can lead to an enhanced phase 
sensitivity. 



A. Mapping to Classical Variables 

The presence of a nonlinearity in the Schrodinger equa- 
tion makes analysis of the dynamics significantly more 
difficult. Notably, the superposition principle is no longer 
valid, making an eigenstate decomposition impossible. 
Indeed, in general there will be more eigenstates than the 
Hilbert space dimension. The system can be approached 
though concepts borrowed from classical mechanics, as 
there is a mapping of our system to a classical Hamil- 
tonian. From this perspective, we study the stationary 
points in the classical phase space, that correspond to 
the eigenstates of the original system. 

Any quantum system with Hilbert space dimension TV 
can be mapped onto a classical system with 2N degrees of 
freedom, namely, the N phases and N amplitudes of the 
state vector. Symmetry under global phase shifts leads to 
a conservation of probability amplitude as an integral of 
motion, and if the Hamiltonian is time-independent, the 
total energy is also an integral of motion. This implies 
that two- mode BECs have non-chaotic dynamics. In the 
three- well case, integrable dynamics are no longer guar- 
anteed. The presence of chaotic dynamics has implica- 
tions for the splitting and recombination elements of the 
interferometer as adiabatic transport through a chaotic 
region of phase space is not possible. Linear quantum 
systems map onto A^-dimensional harmonic oscillators, 
and are necessarily integrable. Introduction of nonlin- 
ear terms leads to more complicated classical dynamics, 
which in some cases can be chaotic. This is indeed the 
case for the three- weh BEC [Ml HQ]. 

We take the classical degrees of freedom to be the 
amplitude squared and phase of the BEC in each well: 
Ni = IV^il^, (/>i = arg (V^i). The number of degrees of free- 
dom can be reduced by two using conservation of proba- 
bility and global phase symmetry: 



N2 = l-Ni-Ns, 

(t>12 = ^1 - ^2, 
^32 = ^3 - ^2- 



(14) 



The classical Hamiltonian is then: 
H/Ml = A(l -Ni- Ns) 

+ ^[Ni+Nl + il-Ni-Nsf 



- 2^12 y/Ni{l -N1-N3) cos (pi2 

- 2^23 VNsil -Ni- N3) cos 032 , 

and the equations of motion are: 

ni2(2iVi+Af3-l)cos(/.i2 . 



(15) 



'12 



y^Ni{l-Ni-N3) 
r223 A^3 cos 032 

ViV3(i - ^1 - ^3) 



g{2Ni + Ns 



Ni/nn = 2^^12x7^1 (1 -Ni- N3) sin 012, 

n23{Ni+2N3 - 1) cos 032 



932 



V^sii - iVi - N2) 

fll2Ni COS 012 

^Ni{l-Ni-^ 



-A 

g{Ni + 2N3 



1), 

(16) 

(17) 



1), 

(18) 

(19) 



Ns/hn = 2i^23\/^^3(l -Ni- Ns) sin 032. 



1. Stationary States 



CTAP, by definition, requires evolution along particu- 



lar eigenstates. In Section [HB] we were able understand 
sensing by the accumulated phase shifting the BEC into 
a superposition of all three eigenstates. However, the 
nonlinearities in the modal BEC approach invalidate the 
superposition principle. Nevertheless studying the states 
of the system still yields useful insight into the structure 
of the resulting phase space. 

The eigenvalues are given by the stationary solutions 
to Eq. (??). To distinguish from the ^ = case we denote 
these stationary states by , with the Dq coinciding with 
the Z^o split state at t = tp/2: 



1 



D',{tp/2) = ^{^, 



^3) = ^ 



50/50- 



(20) 



As mentioned above for the interacting case, the number 
of eigenvalues is no longer limited to three (see Ref. |4T] 
for an in-depth discussion). Fig. [s] shows the eigenvalues 
for various values of ^, using the pulsing sequence defined 
by Eqs. (??) and (??). At higher interaction strengths, 
new eigenstates appear. The appearance of these new 
states marks a bifurcation near the t = {O^tp} limits, 
which disconnects the Dq state from the fully occupied 
states at the end points. Once these eigenstates appear at 
g = gc = A/2± vTTAV4, 50/50 splitting is precluded, 
even in principle. 

The appearance of these additional nonlinear eigen- 
states is a ubiquitous feature of nonlinear systems. How- 
ever, their effect on the dynamics depends on the pulse 
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FIG. 5: Real part of the eigenvalue of the Jacobian at the Dq 
state, (a) A = 0.0. (b) A = 0.1. Transport through regions 
with ^(X) = allows complete fidelity for the protocol. For 
A = 0, the Dq state is unstable for all ^ 7^ 0, albeit with small 



FIG. 3: Eigenenergies for a range of interaction strengths 
with A = 0.1 and tp ^ 00. As the interaction strength is 
increased, new eigenstates appear. When \g\ ^ 1.0, the Dq 
state at t = is no longer ipi, precluding the possibility of a 
50/50 split. 
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FIG. 4: Wavefunction components for the Dq state with non- 
zero ^, A = 0.1 and tp 00. Blue: ipi; green: ip2] red: ips. 
Solid: g = 0.2; long dash: g = —0.6; short dash: g = 1.4; dot- 
dash: g = 2.0. Interactions lead to occupation of the middle 
well during transport, as can be seen for ^ = 0.7 and g = 1.0. 
Note: ipi is real. 



sequence that is employed. In particular Graefe et al [2F 
showed that the extra eigenstates are a permanent fea- 
ture of the eigenspectrum for all gA < and \g\ > |A| 
using a Gaussian pulsing protocol, precluding adiabatic 
transport. Conversely 1^23 / at the start of the proto- 
col, as used here, has the effect of suppressing the emer- 
gence of these extra nonlinear eigenstates for low values 
of 1^1, as seen in the window of robust evolution in Fig.|6] 
This robustness is further considered in the stability anal- 
ysis below. 

From Eq. (??), it possible to calculate the wavefunc- 
tion amplitudes in the three wells as a function of the 
tunnelling rates and nonlinear interaction. Fig. [4] shows 



that for ^ = 0, ?/^2 =0 for all t. However, upon the 
introduction of interactions will in general be non- 
zero. After the appearance of the self-trapped mode, the 
50/50 split state at t/tp = 0.5 is no longer adiabatically 
connected to the fully-occupied well 1 state ai t/tp = 0. 



2. Stability 

We now consider the stability of the Dq state in the 
regime of |^| < \gc\. A stationary point is unstable if 
small perturbations lead to large deviations from the 
point. These deviations make adiabatic transport of the 
system impossible, as keeping the system near Dq re- 
quires that tp ^ 00. 

When 1^1 < l^d, there are just the three eigenstates. 
However, these differ from the noninteracting case, i.e. 
the Dq state's eigenenergy is no longer zero and, in gen- 
eral, '^2(^) 7^ O5 as demonstrated in Fig.|4j The dynamics 
in the neighbourhood of stationary states can be studied 
by making a linear expansion of the equations of motion. 
Solutions are classified by eigenvalues of the Jacobian at 
the stationary point: 



(21) 



where the X'l are the system variables. If the eigenvalues 
are purely imaginary, then the solutions are oscillatory 
around the stationary point. An eigenvalue with a pos- 
itive real part implies the presence of a hyperbolic orbit 
around the eigenstate, which is unstable. 

Fig. [5] shows how the stability of the D'q state changes 
along the pulse sequence for —1 < g < 1, for A = 
(Sja) and A = 0.1 (sfb). For A = 0, the D'q state is 
strictly only stable for ^ = throughout the pulse se- 
quence. However, the introduction of A stabilises the 
D'q state throughout the transport protocol. Specifically, 
for A = 0.1, transport is stable in the range \g\ < 0.2 
and, as can be seen in Fig. [4j complete fidelity for trans- 
port to a split state is possible for extended pulse times. 
At higher |^|, transporting through regions of instability 
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FIG. 6: Fidelity, e, of split as defined in Eq. (??) as a function 
of pulse time, tp, and interaction strength, at t = tp/2 with 
A = 0.1. Regions of good fidelity, e ~ 1, are denoted by 
dark red. The short dashed lines denote g = Qc where the 
bifurcation appears and the long dashed lines denote \g\ — 0.2. 

leads to a loss of fidelity. For an intermediate interac- 
tion strength, this instability is small, and it is possible 
to find a balance between the adiabaticity requirement 
and the slow divergence of the Dq state which still allows 
high fidelity for the splitting protocol. Hence we have 
identified A = 0.1 as advantageous to stable transport. 
We adopt this parameter in the following analysis. 

B. Nonlinear Fidelity of Splitting 

We now explore the robustness of the f-CTAP splitting 
protocol by direct numerical evolution of the nonlinear 
Hamiltonian. Specifically we solve for ipj in the presence 
of nonlinear interactions and finite pulsing times, via: 



We are interested in the transportation of the BEC to 
the split state. We therefore define the fidelity as the 
overlap between the ^50/50 state and the transported 
state, as determined from Eq. (??): 

e= |(M/5o/6ol*(W2))|'- (23) 

Fig. |6] shows the fidelity of BEC splitting in the pres- 
ence of interactions. Qualitatively, for small values of |^|, 
full adiabatic transport is still possible, for tp/2 > 101^"^. 
However, larger interaction strengths lead to oscillations 
and complete loss of fidelity in the splitting procedure. 

The regions of efficient splitting in Fig. [6] can be ex- 
plained in terms of the stability of the Dq state (Fig. [s]). 
For the pulsing scheme given by Eqs. (??) and (??), sta- 
ble transport to the split state is achieved for \g\ < 0.2 
in the adiabatic limit (long dashed lines in Fig. |6|. This 
region is characterised by its close-to-linear behaviour, 
where an increase in total pulse time leads to a corre- 
sponding increase in fidelity. The small instability that 



exists for positive g slightly reduces the fidelity of trans- 
port for long pulse times. The absence of extra nonlinear 
states (Fig. [3| and stability of the Dq state in this range 
(Fig. [5]) means that the fidelity of transport obeys the 
linear adiabatic law up to very long pulse times. 

For 1^1 > 0.2, transport for the Dq state becomes un- 
stable (Fig.|5|, with the growth of this instability increas- 
ing with 1^1 . To achieve good transport in this regime, the 
optimal pulse time must be found. This comes from com- 
petition between the adiabaticity of the transport proto- 
col and the instability timescale of the Dq state, with 
shorter tp/2 being favoured for stronger interactions. At 
9 = Qc-, the appearance of extra nonlinear states near the 
Dq state, which only extend partially along the pulsing 
scheme, as shown in Fig. |3j prevent stable splitting for 
any tp/2. Once these extra states extend to t = tp/2, sta- 
bility in transport is restored. However, the final state 
does not overlap the ^50/50 state. 

The evolution shown in Fig. [6] is the first stage of this 
interferometric process and reversal of the protocol, for 
recombination, will naturally lead to high fidelity trans- 
port analogous to full CTAP transport. Previous work 
by the authors [29 shows that the mapping of CTAP in 
the nonlinear modal approximation to the GPE can be 
achieved with experimentally-realisable trap parameters 
for a three-dimensional BEC containing 2000 atoms. The 
full sensing protocol is the subject of the next section. 

C. Phase Interferometry in the Presence of 
Interactions 

As shown in Section |nB[ the reversal of the f-CTAP 
splitting process gives a phase-sensitive population in 
well 1, Eq. (??). Here we quantify the effect of inter- 
actions on this phase measurement by numerically inte- 
grating Eq. (??) to determine the population of well 1 
at the end of the protocol, Ni{tp -\- r). We assume an 
ideal splitting with Ni(tp/2 + r) = Ns{tp/2 + r) = 0.5, 
£;? = E§ = and = A = 0.1 in Eq. (??), and allow 
the BECs to accrue an initial relative phase difference of 
(j) as in Eq. (??). This investigation of the phase response 
focuses on the behaviour for timescales longer than the 
linear adiabatic pulse time, where tp > 100Q~^. 

The final population of well 1 for tp = 200Q~^^ tp = 
bOOn-^, tp = 1000(^-1 and tp = 20001^"^ (Fig.[7|) exhibit 
large variations that prevent a robust phase measurement 
for 1^1 > l^cl- For \g\ < 0.2, the phase response is similar 
to the noninteracting limit [Eq. (??)] and coincides with 
the region of stability for the Dq state (Fig. |5|. 

We have computed first-order corrections to the peri- 
odic orbits using classical perturbation theory For 
small 1^1, the deviation of the final population well 1 from 
the noninteracting case, AA^i, is given by: 

(.4) 

which agrees with the numerical simulations in Fig. [S] 
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FIG. 7: Occupation of well 1 after ideal BEG splitting and recombination using reversal of f-GTAP protocol for (a) tp — 200^1"^, 
(b) tp — 500r2~^, (c) tp — 100011"^, and (d) tp — 2000^"^ with A = 0.1, as a function of the phase difference (0) between 
wells 1 and 3 at t — tp/2 and the strength of the inter- atomic interactions \g\ < 1. 
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FIG. 8: Deviation in the population of well 1 from the noninteracting expected value, A^id^^l < 0.1) — Ni{g = 0), as a 
function of the phase difference (0) between wells 1 and 3 at t = tp/2, for (a) tp = lOOQ"^, and (b) tp = 2000^"^, where 
Ni{g = 0) = cos^((/)/2). The classical first order perturbation, Eq. (??), is plotted for \g\ — 0.05: dashed curves, and \g\ — 0.1: 
solid black curves. 



The non-zero /S.Ni implies an increased phase sensitivity 
due to the interactions, which is robust with respect to 
the pulse time. 

Consistent with Eq. (??), we numerically find that 
AA/"! is zero for cj) = {0,7r} (Fig. |8|. For long pulse 
times, Fig.jsjb), the longer time spent traversing the un- 



stable regions of phase space leads to a larger deviation 
near (j) = {0,7r}. For small attractive interactions. Fig. [s] 
shows that for ^ the interferometer is stable. Con- 
versely, for small repulsive interactions. Fig. [8] shows that 
for (/) ^ TT the interferometer is stable. 

Even though bifurcations do not emerge till \g\ > 
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\gc\ ^ 1.0, Fig. [t] demonstrates that for interaction 
strengths \g\ > 0.5, the instabihty of the symmetric and 
antisymmetric stationary states at (/) = {0,7r} lead to 
large variations in density. Large density fluctuations re- 
duce the possibility of recovering any phase information 
across all of ^ for ^ > 0.7. 

In Fig. we show the regions where phase-sensitive 
measurements can be performed in the presence of non- 
linear interactions, i.e. \g\ < 0.35. Since nonlinear 
atom-atom interactions can introduce unwanted density- 
dependent fluctuations, which compromise phase sensi- 
tivity, atomic interferometry experiments often reduce 
these interactions via a Feshbach resonance [43l, '44| . 

In the small interaction limit with \g\ < 0.1 a robust 
phase-sensitive measurement of Ni can be made by ac- 
counting for the correction given by Eq. (??). In this 
regime, Fig. |9] (dark blue lines) shows a smooth response 
of Ni with respect to phase, with deviations for the non- 
interacting result shown in Fig. [Sj The instability of the 
(j) = 7T state for attractive interactions, seen as deviation 
from the first-order correction in Fig. [S) does decrease 
the pulse time sensitivity relative to the stable regions. 
However, the f-CTAP interferometer should be robust for 
a broader range of interactions. 



1. Enhanced phase sensing 

Here we highlight how a f-CTAP based interferometer 
has two distinct regimes of enhanced phase-sensitive op- 
eration, facilitated by control of the interaction strength 
via a Feshbach resonance. Regime 1 is characterised by 
sensing around ^ = {0, tt} for short pulse times whereas 
regime 2 requires much longer pulse times, but allows for 
tuneable sensing. 

In regime 1: For shorter pulse times [see Fig. [7]^a) and 
^h)] around (j) = {0,7r}, an enhanced phase-sensitive 
measurement of A^i can be made. For small |^|, the 
first-order corrections to the final population in well 1 
[Eq. (??)] are small, and there is very little sensing en- 
hancement. However, for greater |^|, large but stable 
oscillations appear in these regions, leading to enhanced 
sensing. This occurs for 0.2 < g < 0.6 {(j) ^ 0) and 
-0.2 > g > -0.6 {(j) ^ tt) [Fig. ^d^)]. The emer- 
gence of these oscillations can be seen in Fig. [8]^a) where 
= {0, tt} oscillations start to occur for ^ < and ^ > 
respectively. As can be seen in Figs. [7|^a)-[7]^d) the range 
of interaction strengths where the (j) = {0, tt} states re- 
turn to A^i = {1,0} reduces with increasing pulse times. 

In regime 2: A highly enhanced phase-sensitive mea- 
surement of A^i can be made at the point indicated by 
the dashed line in Figs, [ojc) and[9]^d). This stems from 
the large gradient in the population phase response in the 
vicinity of the dashed lines. The position of this bound- 
ary could be controlled by tuning the interaction strength 
via a Feshbach resonance. This would, for instance, allow 
the position of maximum sensitivity for a given phase to 
be changed depending on the application. This opens up 



the possibility of using a feedback-based phase measure- 
ment to ensure maximum sensitivity to particular phase 
changes. 



2. Robustness 



As shown in Section |III A2[ the adiabaticity criteria, 
Eq. (??), is invalidated by the introduction of the non- 
linear interaction term. Hence, here we associate robust 
phase-sensitive measurement with a high sensitivity to 
changes in relative phase, dN/dcj)^ and a low sensitivity 
to changes in pulse time, dN/dtp. 

For the small interaction limit, where the di- 
vergence due to the instability is relatively slow. 
Figs. [lO]^a) and [lO|b) shows that as we increase the 
pulse time from tp = 100^1" ^ to tp = 2000^1" ^ the pulse 
time sensitivity decreases by a factor of 10~^. However, 
this deviation is small with respect to phase sensitiv- 
ity. While this regime exhibits minimal enhancement of 
phase-sensitivity, it will have the robustness and stability 
of a noninteracting f-CTAP interferometer. 

In regime 1: there is some reduction in pulse-time sen- 



sitivity relative to the small interaction limit [Fig. 10 



however the pulse time sensitivity is still a factor of 10 less 
than the phase sensitivity for tp = 100^1"^ and is further 
reduced by a factor of 10^ for tp = 20001^-^ [Figs. [To|;c) 



and 10 ^d)]. The nonlinear interaction also leads to a 
smaller enhancement of the phase-measurement for g < 
0.3 and g > 0.3 for <p = {0,7r}, which is more robust to 
changes in pulse time due to the reduced instability. De- 
spite this. Fig. 10 ^d) shows that the phase enhancement 



persists for longer pulse times. 

In regime 2: Surprisingly, the robustness of sensing 
for this range of phase shifts is increased for longer pulse 
times where peaks in phase sensitivity no longer align 
with peaks in pulse time sensitivity [Fig. fTo|d)]. 



IV. CONCLUSION 

We have demonstrated that the f-CTAP protocol is 
a valid methodology for the coherent spatial splitting 
of a BEC in the presence of interactions. Reversal of 
the f-CTAP protocol provides a robust phase sensitive 
measurement as an alternative to traditional methods. 
This interferometric process is robust to changes in pulse 
time once in the adiabatic regime for the small interac- 
tion limit. Also, even though the adiabatic principle is 
not strictly valid in this nonlinear system, it is possi- 
ble to find splitting and recombination timescales that 
achieve the balance between maintaining the system in 
the Dq stationary state and divergence due to instability. 
Phase-sensitive measurements in the presence of nonlin- 
ear atomic interactions can lead to an enhancement in 
sensitivity without significant loss in robustness with re- 
spect to changes in pulse time. 




FIG. 9: Phase sensitive measurement of the occupation of weh 1, Ni{tp)/NT, for interatomic interaction strengths \g\ < 0.35, 
after recombination times of tp — 100Q~^, (a) and (b), and tp — 3000O~^ , (c) and (d). Ideal sphtting of a BEG as defined in 
Eq. (??) is assumed. 
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